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1. Introduction 

According to Scliafke and Sclimidt P the generalized Heun equation (GHE) is a second 
order differential equation of the form 

y{z) + [ + r + - + a]y{z) + — — —y{z) = 1.1 

\ z z—1 z — a J zyz — \)yz — a) 

where a G C\{0, 1} and yUo, /^i, [i'2-, /^O; A, are arbitrary complex numbers. 
Moreover, 0, 1 and a are simple singularities with exponents {0, /io}, {0, /^i} and {0, yU2} 
respectively, while oo is at most an irregular singularity of rank 1. Following pp we 
choose parameters A := (Ai,A2,A3) G such that the last coefficient of ()1.1|) takes 
the form 



;3o + ^iz + ^2Z^ _ 1 fl- fiA f I- fip\ , (a/2)(l - /ifc) + A 



with 2:9 := 0, Zi := 1, ^2 := 0^ G C\{0, 1}. Note that for fixed /zq, /ii, /i2, « and a A is 
uniquely determined by {Po, Pi, P2) and vice versa. 

To underline the importance of equation (jl.lj) we recall that it contains the ellipsoidal 
wave equation as well as the Heun equation {a = P2 = 0) and thus the Mathieu, 
spheroidal. Lame, Whittaker-Hill and Ince equations as special cases. The aim of our 
work is to present some specific examples where the GHE arises in physics with particular 
attention to QFT theory in curved space-times. 

Schafke and Schmidt obtained under certain conditions on a all connection coefficients 
between the Floquet solutions at the singularities 0, 1 and a, thus determining the full 
monodromy group of (jl.lll . Furthermore, they showed that the connection coefficients 
can be computed by four-term recurrence relations. 

Notice that for a = and /32 = Ai + A2 + A3 = the equation (jLip becomes the Heun 
equation where the point at infinity is a simple singularity with exponents 

+ (l/2)(/io + /ii + At2), -z/ + 1 - (l/2)(/io + /ii + /X2)} 

and 
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7 (/io + /i? + /i2) - ^ - - aA2 



2. The Dirac equation in the Kerr-Newman metric and the GHE 

In Boyer-Lindquist coordinates (t, r, t?, with r>0, 0<'i9<7r, 0<(/9<27r the 
Kerr-Newman metric is given by j2] 

,a A-a^sin^?? ,n 2a sin^ t^fr^ + - A) ,,,.,,2,, 
ds^ = dt^ + ^ dt dip - (r^ + a^y sin^ i9- dcp^ 

- ^ dr=^ - S d^^ 
A 

with 

S := S(r, e)=r'^ + cos^ 9, A := A(r) = - 2Mr + + 
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and 



E := E(r, ^9) = 1 - a'Y{r) sin^ 1^ 



7(r) 



+ 



where M, a and Q are the mass, the angular momentum per unit mass and the charge 
of the black hole, respectively. In the non-extreme case > + the function A 
has two distinct zeros, namely, 

r_ = M - ^/M^-a 



r+ = M+ 



the first one corresponding to the Cauchy horizon and the second one to the event 
horizon. Moreover, A > for r > r+. In the extreme case = + the Cauchy 
horizon and the event horizon coincide since A has a double root at ri^ = M. According 
to Penrose and Rindler [3] the Dirac equation coupled to a general gravitational field V 
is given in terms of two-component spinors (0"^, ) by 



VTLf. 



{V\ - leV^ )xa' 



where we used Planck units h = c = G = 1. Furthermore, V^^' is the symbol for 
covariant differentiation, e is the charge or coupling constant of the Dirac particle to 
the vector field V and rrie is the particle mass. The Dirac equation in the Kerr-Newman 
geometry was computed and separated by Page |3] with the help of the Kinnersley tetrad 
[H]. In view of the separation of the Dirac equation we choose to work with the Carter 
tetrad . If we make the ansatz [71 IH] 



iXi' 



R+{r)S4^) 



k = k + 



1 



where u and G Z denote the energy and the azimuthal quantum number of the particle, 
respectively, and 5* is the non singular matrix 

S{r,d) = A"4 diag 



^/r — ia cos 'd \/r — iacos-i? -y/r + ia cos A/r + ia co^'d , 
then the Dirac equation decouples into the following systems of linear first order 
differential equations for the radial R± and angular components S± of the spinor \l/ 




^AI?_ 
\mf.r — A 

—L- X + arUe cos 6 
X — anip cos 9 C+ 



0, 




where 



d .K{r) 
d^ 



K{r) = u}{r'^ + a?) — eQr + ak, 



£± = + ^cot^±Q{^), Q{^) = ausin^ + kcsc^. 

CI t/ ^ 



(2.1) 
(2.2) 

(2.3) 
(2.4) 
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For a detailed derivation of (|2.ip and (|2.2j) we refer to [H]. Let us briefly recall that 
in the case a = the components S±. of the angular eigenf unctions can be expressed 
in terms of spin-weighted spherical harmonics [HI El whereas for a 7^ they satisfy a 
generalized Heun equation JT]. We show now how to transform the radial equations 
for the components R± of the radial functions in such a way that they are reduced to a 
GHE. In what follows a prime ' always denotes differentiation with respect to r. First 
of all let us bring ()2.3|) into the form 

i?'_(r) + z/(r)i?_(r) = f{r)R+{r) (2.5) 
Fljr) + V{r)R+{r) = J{r)R4r) (2.6) 



with 



A(r) ' ^/AW ■ 

From the above system and the definitions for u and / it can be immediately verified 
that R± = R^. As a consequence we just need to derive the ODE satisfied by A 
simple computation gives 

i?" + pR'_ + qR- = 0, RI+ pR!^ + qR+ = 

with 

where to keep notation simple we write R±, p, q, f, v instead of R±{r) and so on. 
Notice that g — cu^ — ml for r +00 as can be seen easily from the expansion in 
partial fractions below. Moreover, we shall restrict our attention to the case |c<j| > mg. 
For future convenience we write p and q in terms of partial fractions as follows 

Po , P- , P+ 



p{r) 



where 



and 



where 



1 1 . A 

Po = -1, p-=p+ = -, ro = 1 — 

2 rrif. 



, , „ Tn XL XV 5'n S— S-L 

q(r) = X+^^ + + ^^+ '-^ + -, ^ + 7 ^ 

r — Tq r — r_ r — [r — ro)^ [r — r_)"^ [r — ? 



X 
To 

X_ 



2 2 



iA 



iB 



r+ - ro r+ - ro 

2AB - my_ - iuj i{A + B) 
r>—r^ 2 2(r I — r_ 



2Buj' 



iB 

) r_ - ro ' 
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2 2AB-mlrl-\'^ iu i{A + B) iA 



2Auj' + ^ + 



r+ — r_ 2 2(r+— r_) r+ — tq' 



^0 =0, 5_ = 5^ - ^, ^+ = - y 



with 



Finally, by transforming R_ according to 

R_{r) = {r - r^y- {r - r+y+e'^^'R^ir) 
and introducing the new variable z G C defined by 



r I — r_ ' 



it can be checked that if l± and 7 are such that 
7 = ±1^/^0-2 - ml, 

4 + ip± - 1)^± + s± = o, 

then (ESI) reduces to jHH) with 

11 , 
/io = --2i_, /ii = --2i+, /i2 = 2, a = 2(r+-r_)7, 

Ao = (r+ - r_)T_, Ai = (r+ - r_)T+, A2 = (r+ - r„)Tc 







1 M X/rrie 
7: + + i- 



r_|_ — r_ — T- 

3. Static perturbations of the non extremal Reissner-Nordstrom solution 
and the GHE 

Let us consider the following Lagrangian describing four-dimensional gravity coupled 
to a U{1) gauge field and a real massive scalar, with a non-renormalizable but gauge 
invariant coupling between the gauge field and the scalar, namely 



2 - ' " 1 + 

where R and Gat are the scalar curvature and the Newton constant, respectively, and i 
is a parameter having the dimensions of length. In order to construct static solutions 
with magnetic charge g we require that 

ds'^ = gudf + grrdr'^ + r^(d6'^ + sin^ 6 dip'^), 
1 

2' 

= Hr)- 



F2 = -F^^ dx^ A dx" = A-ng = g d9 AsinO d(p. 



Moreover 



By setting M« = 1/V8^FG^ = 1 and 

we obtain for and i? the following system of equations [12] 

," /2 ^, 1 ,'2\ oFidVeff 

+ - - 25 + -r0 » - 
\ r 2 



2B 



1 ,/2 r,R^ ^ 2i? 1 — e 

-0 + e^^K// - — + = 0. 

Linearizing the above equations of motion around = and requiring that the horizon 
is at th = 1 yields 



+ (?-2B'V'=e-^(0.rM, 



r 



B{r) = --\ogil-— + ^], -—^{0,r)=m^-^ 



2 ^ V 47rr 2r2 / ' 902 ^ ' ' 



ocl + ^^^ ^(r-l) for r^l, 



oc for r cxD 

r 

where M is the mass of the black hole. We refer to ^2] for the determination of M by 
fitting e~2'^ to the Reissner- Nordstrom form for large r. Taking into account that the 
requirement th = ^ is equivalent to require that M = An + 2g'^ the ODE for can be 
brought into the form 

+p(r)0 + g(r)0 = 

with 

( . Ph ^ P~ 2 ^2 

P[n = H , Ph = , P- = -r-PH, r_ = — , 

r — 1 r — r_ 1 — r_ 2tc 

g(r) = -m + ^ H ^ H , 

9 .2 A ^ M ^^^^^ - ( 2 ?00 

goo = 27r£ , go = H ^oo, 9h = , g_ = m r. 



Finally, by means of the transformation 

0(r) = r"e-'""V^(r) 

with 5 such that 52 — 5 + goo = it can be checked that ip satisfies the GHE (ll.lj) with 
/io = 1 — 25, /ii = 1 — /^2 = 1 — « = —2m, 
Ao = go, Ai = g^, A2 = g-. 
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4. Conclusions 

In this paper we have related the radial Dirac equation for a massive fermion in the 
Kerr-Newman geometry to a GHE. In addition, we found that the GHE describes also 
static perturbations for the non-extremal Reissner-Nordstrom solution. We reserve the 
study of other examples where the GHE appears for future investigations. 
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